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Abstract 

Measures generated by Iterated Function Systems composed of un- 
countably many one-dimensional afFine maps are studied. We present 



I— I • numerical techniques as well as rigorous results that establish whether 

pH ' these measures are absolutely or singular continuous. 
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> ' 1 Introduction 

in 

O^ ■ In this paper I want to describe an example of the fruitful interplay between 

1_" \ the theory of dynamical systems and numerical analysis: I want to show how 



theoretical questions such as singularity (or continuity) with respect to Lebesgue 
of a dynamical measure can be attacked from a numerical point of view, and 
viceversa how particular features in the numerical analysis of measures can be 
properly explained by concepts in the theory of dynamical systems. 

Let us consider the iteration of maps 4>\ : X — > X from a compact metric 

space X to itself, labelled by the variable A which belongs to a measure space A, 

j^ ■ on which the probability measure a is given. This process gives rise to what is 

called an Iterated Function System, or I.F.S. , [101 [Tj with invariant measure /i, 
which can be defined as follows. Consider the transfer operator T on the space 
C(X) of continuous functions on X, via 

(Th)(x)= [da(\)(ho<f> x )(x). (1) 

Then, let T* be the adjoint operator in the space of regular Borel measures on 
X. An invariant measure of the I.F.S. is the fixed point of T*, 

T*(») = n. (2) 
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Suitable hypotheses can be formulated in order for \i to be unique, given the 
choice of the set of maps <j)\ and of the distribution a [18] . We are interested in 
the nature of /z: is it of pure type? If so, is it pure point, singular, or absolutely 
continuous? What characteristics of cj)\ and a have importance in this regard? 
This problem belongs to a classical topic of research in dynamical systems, that 
looks for a.c.i.m., that is, absolutely continuous invariant measures. Indeed, 
it is my opinion that singular continuous measures are equally, if not more 
interesting, in many respects. 

The plan of this paper is to study this problem in a class of maps of the real 
line. We will derive algorithms to reduce the quest for an invariant measure to a 
fixed point problem. If in addition there exists an attractive fixed point in a suit- 
able space of densities, absolute continuity of the measure will follow. We will 
also describe the Fourier transform of the measure /i, its Mellin transform and 
its Sobolev dimension, that will also lead to numerical and theoretical methods 
to determine absolute continuity or singularity of the invariant measure. 

2 Infinite Affine Iterated Function Systems 

Iterated Function Systems have been originally introduced and studied in [TUlfT]. 
although in some form their history goes much further back in time, see |21j . 
They have become a versatile mathematical tool with applications to image 
compression [3] [T^] , quantum dynamics [TS] and much more [5] . In its simplest 
form, an I.F.S. is a finite collection of contractive maps of a space X into itself: 
when iterated randomly, these maps produce a stochastic process in X with 
invariant measure /i. Interesting results and applications are found for affine 
one dimensional maps of the kind 

Mx)=5{x-P)+p. (3) 

Each of these maps has the fixed point j3 and contractions ratio < 6 < 1. 
When dealing with finite collections of affine one-dimensional maps, the problem 
of constructing their Jacobi matrix, and therefore of computing integrals with 
respect to these measures, has been solved in [T31I71IT5]. 

Rather than working with a finite set of maps associated with the pairs 
(S, j3) , as done usually, in this paper we adopt a generalization by which S is 
fixed to a positive constant value, strictly less than one, while j3 is free to vary 
into a finite interval of the real line, that, without loss of generality, we may 
understand to be [—1, 1]. This choice is called an infinite, homogeneous I.F.S.. 
It has been originally investigated by Elton and Yan [Bj. 

Using the maps 4>s,p we can define a stochastic process in X = [—1,1] via 
the following rule: given an initial point x G [—1,1], choose a value of j3 at 
random in [— 1, 1], according to a distribution <r(/3) (whose support may contain 
an infinity of points) and apply <f>s,f3 to map x into (frs,p(x)- Iterate the procedure. 
A general theorem, due to Mendivil [18], guarantees that there exists a unique 
invariant measure \i for this stochastic process. In addition, this measure can be 
found, probability one, by the Cesaro average of atomic measures at the points 



Xj of a trajectory of the process: - X^j ~~ ^ M- I n a rather pictorial view, 
we may describe the process above by saying that the point x is the location 
of a predator in a rather peculiar chase of the prey, located at the point j3: 
the predator moves towards the prey, but as soon as its distance from /3 is 
reduced by a factor 5 the prey disappears, to reappear instantaneously at a new 
location, and the process is repeated. Therefore, a is the distribution of the 
position of the prey, that, along with the value of S yields the distribution ^ 
of the predator's positions. The measure \i can be equivalently defined by eq. 
((2]): its approximation properties and the related inverse problems have been 
discussed in [3j Q2] . 

The transfer operator T for infinite, affine homogeneous I.F.S. takes the 
following form: 

(Tf)(x)= f da(/3) f o ^(x) = [da((3)f(Sx + Sp), (4) 



where / is any continuous function, and where, to simplify the notation, we 
have introduced the symbol 6 := 1 — 6, to be used throughout the paper. Eq. 
12) now becomes: 



d(T*n)(x)f(x) := J d»(x)(Tf)(x) = j dfi(x) j da(/3)f(Sx + 6/3). (5) 
This is equivalent to say that, if fi is invariant, then the equality 

dn(x)f(x)= [ dii(z)(Tf)(x) (6) 



holds for any continuous function /. In this setting, our problem can therefore be 
formulated by asking what characteristics of the constant 5 and of the measure 
a influence the nature of (i. We start attacking this question in the next section, 
in the case when a is the Lebesgue measure. 

3 Test Case I: a is the Lebesgue Measure on 

[-14] 

In this section, we consider a case that can be treated analytically as well as 
numerically. It is defined by letting a be the Lebesque measure on [—1, 1]. By 
proving theoretically the right answer for the spectral type of fj,, we shall be 
able to use it as a test of the numerical techniques that we shall introduce in 
the sequel. We first prove a general result 

Lemma 3.1 Suppose that a is absolutely continuous with a bounded density. 
Then, so is fj,. 

Proof. Apply the balance relation (J6j) to f(x) — Xb c (s){ x ), the characteristic 
function of B e (s), the ball of radius e centered at s (technically, eq. (J2]) holds 



also for summablc functions), to get: 

p(B e (s)) = J dp(x)a(B e/S ( S —J^)). (7) 

If a has a bounded density, this means that the following limit exists and is 
equal to p(cx; y), the density of the measure a at the point y: 

HmJ-a(B e {u))=p(* ]y ). (8) 

Using this dehnition, we can now divide both sides of eq. (O by 2e and take the 
limit for e — > 0. Thanks to the dominated convergence theorem the limit can 
be taken inside the integral at r.h.s. to obtain 

If s — ox 

P(W a) = j / dp(x)p(a; — — ). (9) 

From this last relation the thesis easily follows. I 

We now specialize the theory to the case when a is the Lebesgue measure, 
that is to say, all possible values of j3 in [—1, 1] are equally probable. 

Theorem 1 Let < 8 < 1 and let a be given by p(a;j3) — \x{-i,i\{P)- Then, 
the invariant measure p is absolutely continuous, with a density p(p; x) that is 
infinitely differentiable and non-analytic. 

Proof. The fact that p has a bounded density follows from Lemma 13.11 Fur- 
thermore, inserting the value of p(a; [3) into eq. (|9]), we obtain 

p(p;x) = ±p([^,^-}), (10) 

and, taking the derivative with respect to x, 

p'(^x) = —^[p(p;—^—)-p(p;—^—)]. (11) 

This equation can be iterated to show the existence of derivatives of p{p; x) of 
all orders. Furthermore, observe that, when taking x = — 1, both points ( ~g~ ) 
and ( -1 f ) are to the left of —1, that is, outside of the support of p, so that all 
derivatives of p(p;x) are null in x = — 1. Therefore, p(p;x) cannot be analytic. 
I 

In the next section, we will describe a numerical technique to compute the 
density of this invariant measure. 

4 A Density Mapping 

Let us suppose that the invariant measure p of an infinite, affine, homogeneous 
I.F.S. is absolutely continuous, with density p(p;x). We want to derive a nu- 
merical technique to compute this density. In the case of Bernoulli I.F.S., to be 



discussed later in section [5] this idea has been proposed in [3] , without quanti- 
tative tests of convergence. Observe for starters that the action of the adjoint 
operator T* can be transferred on densities. Use f(x) — Xb c (s)( x ) m ec L- © to 
get, for any probability measure v. 

(T*i/)(B.(aO) = J da((3)v(B e/s (^l)). (12) 

Clearly, a result quite similar to Lemma 13. II holds: 

Lemma 4.1 Suppose that v is absolutely continuous with a bounded density. 
Then, so is T*v . 



In fact, operating as in the proof of Lemma 13. II of the previous section, one also 
finds the density of TV: 

p(T*v; x) = \J da(l3)p(v, ^^). (13) 

Eq. (|13[) is particularly suited for analytical computations. It can also be 
implemented numerically, whenever a convenient representation of p is found, 
and integration with respect to a is numerically feasible. In this perspective, it 
is best to consider a Fourier space representation. Take therefore f(x) = e~ wx 
in eq. (O and use the notation 



u(y) := / dv{x)e~ tyx (14) 

to indicate the Fourier transform of an arbitrary measure v, to get 

f%{y) = p(5y) ■ a(Sy). (15) 

This is a rather crucial relation, already discussed in [6^, that links the Fourier 
transforms of a and /i. We shall make use of it repeatedly. Observe also that 
when v is absolutely continuous, that is, dv{x) = p(v;x)dx, i>{y) in eq. (|14j) can 
equally well be seen as the Fourier transform of the density p(y; x). Furthermore, 
since the support of p, is enclosed in [—1,1], we consider the Fourier coefficients 
of p(p;x) over the basis set e~ l7Tkx , with k integer. Let us call them Ck(p)- 

c k (p) :=/2(7rfc) := [ dp(x)e- lnkx = [ e~ l7Tkx p(p; x)dx. (16) 



Setting y = irk yields the Fourier coefficients of T* p at l.h.s. of eq. ([15]). Simple 
formal manipulations of the right hand side prove the Lemma: 

Lemma 4.2 The Fourier coefficients of T* p depend linearly on those of p via 
the following operator 

oo 

Cfc (T*p) — a(Snk) \_, sinc(n(j — kS)) Cj(p). (17) 

j=— oo 



This lemma is theoretically inspiring, and, at the same time, it can be turned 
into a computational procedure. From the theoretical side, consider the fol- 
lowing argument. Suppose that we start from an initial distribution /iQ whose 
Fourier coefficients decay extremely fast for |fc| larger than some value, say ko- 
As it is well known, this means that (1q is absolutely continuous, with a very 
regular density. We observe from eq. (TTT1) that this distribution of coefficients 
may spread under the action of T* . In fact, two phenomena are competing, 
in this regard: On the one hand, Cfc(T*/uo) receives a contribution from Cj(/j,o) 
that is maximum for j ~ kS, and therefore tends to "populate" higher frequen- 
cies, that is, to decrease the regularity of the densities. On the other hand, 
this contribution is multiplied by a(Snk), which tends to zero as k grows, if a 
is sufficiently regular, and therefore its effect is diminished. On this basis, a 
theoretical estimate can be carried out, in line with the project described in the 
introduction: use numerical analysis to inspire theoretical proofs. This estimate 
will be presented elsewhere. Instead, we now proceed with numerical techniques 
and experimentations, that illustrate this point. 

We use eq. (fTT|) as in a fixed point method: that is, we define a sequence 
of measures \i n := T* n /io, with /j,q an arbitrary starting measure and we iterate 
eq. (|I7|). This leads to 

Algorithm 1. 

Fix an integer size M and a threshold r\ and compute the vector of Fourier 
coefficients a(5wk), for k = — M, . . . , M. 

1 Initialization: Put n — 0. Choose a suitable density to define the initial 
measure no and its Fourier coefficients Ck(no), for for k = —M, . . . , M. 

2 Iteration: given the Fourier coefficients of /i n , compute Ck{T*fi n ) = Ck{fi n +i) 
via eq. ([FT)) , restricting the summation over j to the range [— M, M]. 

3 Control and termination. Compute a distance (chosen among various 
possibilities, see later) between fj, n and p> n +i- If difference is less than the 
threshold 77, stop. Otherwise, augment n to n + I and loop back to 1. 

In numerical experimentations we use as initial measure fj,o either the uniform 
density of [— 1, 1], given by p(no;x) = ^X[-i,i]( x )> whose Fourier transform is 
readily computed, or the Fourier transform of a Gaussian distribution centered 
at zero, of variance v2/w, with uj < 1. 

To test this technique we apply it first to the case where a is the Lebesgue 
measure on [—1,1], discussed in section[3l for which we have a precise theoretical 
result. In Fig. Q] we plot the density of the measures obtained by the first five 
iterations of the method, p(fi n ; x), n — 1, . . . , 5, when and 5=1 and M = 100. 
The last two curves are practically undistinguishable, indicating convergence. 
Indeed, a distance function is required to assess this fact with rigour. Various 
choices are possible. In Fig. [2]we plot all of the following distances, in the same 




Figure 1: Densities p(p n ;x) in the iteration of Algorithm 1, for n = 1, . . 
when a is the Lebesgue measure. The initial density p(po;x) is constant. 
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case of Fig. [TJ but on a larger number of iterations. 
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In all cases we observe exponential decrease of di{p(ji n ;x),p(jj, n —i;x)) with n, 
which implies convergence to a density pm(p~,x) that is an approximation of 
p(p;x). Validity of the approximation can be checked by increasing the value 
of M, see the next Section and Sect. EJ where a drastically different case is 
encountered. 



5 Test Case II: a is a two-atoms Bernoulli Mea- 
sure 

Opposite to the previous choice of the measure a, lies the Bernoulli measure, 
that is, an atomic measure composed of two atoms. This is the first non-trivial 
case, because one atom alone leads to the equality p = a. Recall the predator- 
prey interpretation: if the prey at /3 does not move, the predator x converges 
to it, and the invariant distribution p is an atomic measure located at the same 
point. Let therefore a be defined as the sum of two atomic measures located at 
1 and -1: a{(3) = \{^M + A^/3)), or / f(p)da{0) = |(/(1) + /(-l)), for 
any continuous /. 

The prey appears with equal probability (this symmetry can be easily changed, 
though) at either end of the interval. Contrary to what it might seem, this prob- 
lem is far from being trivial. It has a rich and classical history, going back at 
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Figure 2: Distances di(p(p n ; x), p(p n —i; x)) between successive densities in the 
iteration of Algorithm 1. Symbols for I = 4 are red lines, crosses; I = 3 blue 
lines, stars; I = 2 green lines, stars; and I — 1 magenta lines, squares. 



least to Wiener, Wintner, Erdos and others in the 1930's. Present knowledge 
can be summed up in the 

Theorem 2 Let c(/3) = |(A_i(/3) + Ai (/?)). For any 5, p is of pure type. It 
is singular continuous, supported on a Cantor set, if 5 < i. When 6 = \, p is 
the Lebesque measure on [—1,1]. There exist two constructive, countable sets of 
values of 6 > i for which p is singular continuous, and absolutely continuous, 
respectively. For almost all S > \, p is absolutely continuous. 

It is just the case to remark here that this theorem collect significant results 
obtained along more than seventy years of research, see [H Q]5] for lists of ref- 
erences. It is easy to understand what happens in the case S < \. Consider 



cq. ([T3f and let p(po) be the uniform density on [— 1, 1]. It follows by direct 
computation that p{p n ) is a piece-wise constant function for any n, which takes 
values 2~ n ~ 1 5~ n on a set of 2™ disjoint intervals of equal length 28 n , and zero 
otherwise. These intervals constitute the usual generations in the hierarchical 
construction of a Cantor set (take 6 = | to obtain the classical, ternary Cantor 
set). Clearly, the sequence p(p n ) does not tend to any function and this is a 
remarkable example where measures converge while densities do not. 

The bounded variation norm is particularly suited to illustrate this point: 
we have that \\p(p n )\\ B v = 5~ n and \\p(p n ) - p{p n -i)\\BV = 26~™(^j - 1). 
Both quantities diverge as n — > oo. This provides a second test for algorithm 
I, that should not converge in this case. In numerical experiments we indeed 
observe an exponential increase in the bounded variation norm ||p(^n)||sy, that 
only saturates because of the finite cardinality 2AI + 1 of the basis set. Quite 
obviously, increasing M pushes the saturation point to the right: higher and 
higher frequencies are required to describe the densities p(p n )- This much for 



Algorithm I: a better technique will be described momentarily. 

We need therefore to face numerically the possibility of fi not being abso- 
lutely continuous, with slowly decaying (or not decaying at all!) Fourier coeffi- 
cients. In the next section we tackle this problem, first theoretically and then 
numerically. 

6 Fourier Transforms and Sobolev Dimension 

Let us therefore develop techniques for the case when \x may be singular contin- 
uous. A few general results must be quoted at this point. Consider the Mellin 
transform of a function h defined on [l,oo]: 

/oo 
V^Ky) dy. (19) 

The integral may diverge if Re(z) is too large. The supremum of the set of 
values of Re{z) for which one has convergence is called the divergence abscissa 
of the Mellin transform. Define the Sobolev dimension of a measure \i as the 
divergence abscissa of the Mellin transform of h(y) — |/2(y)| 2 : 

d s (fi) := sup{s e R s.t. Mi(|/i(-)| 2 ; s) < oo}. (20) 

Clearly, d B (jx) > 0. When d s (fJ,) is less than, or equal to one, it coincides 
with £>2 (/•*)> the correlation dimension of the measure fi, a common quantity in 
the multifractal analysis of measures 20 . In this case, the leading asymptotic 
behavior for large y of the Cesaro average of h(y) = \fi{y)\ 2 is y~ D2 , see QUE] 
for a discussion of this and other asymptotic behaviors of singular continuous 
measures. To the contrary, when d s (fi) > 1, one knows that fj, is absolutely 
continuous and its density p(/j,) has fractional derivative of order (d — l)/2 in 
L 2 (R) for all 1 < d < d s (/x). If in addition d s (fi) > 2, then p(/j,; x) is a continuous 
function, the larger d s the higher its regularity. To sum up, the regularity of a 
measure can be assessed by the study of the asymptotic behavior of its Fourier 
transform. 

To do this in our case, observe that equation (|T5|> can be iterated, to get 

n— l 

£(t/) := T"So(2/) = MS n y) n fifth)- ( 21 ) 

j=o 

Observe furthermore that, as n tends to infinity, flo(S n y) tends to /io(0) = 1 for 
any y, so that this proves the 

Lemma 6.1 The Fourier transform of fi n can be computed according to eq. 
(21]) , and that of fi in the form of the infinite product 

oo 

m = l[&(Fh)- (22) 

3=0 



Therefore, y is an infinite convolution product of rescaled copies of the measure 
a, a fact well known when a is the Bernoulli measure and fl(y) is an infinite 
product of trigonometric functions. This observation is the basis of many the- 
oretical investigations. It also entails a numerical technique to compute either 
JMi(y) (just use eq. (glj) or p,(y): 

Algorithm 2. 

Fix a threshold r\ > and y £ R. 

1 Initialization: Put n = and let (f>Q = 1. 

2 Iteration: compute %p n = a(S n Sy). Update 4> n +\ = (finrfn 

3 Control and termination. Compute \ip n — 1|. If this difference is less than 
the threshold rj, stop. Otherwise, augment n to n + 1 and loop back to 2. 

At this point, we can outline a two-fold strategy for detecting numerically 
the continuity properties of a measure p generated by affine I.F.S.. In the first 
approach, we use Algorithm 2 to compute the Fourier coefficients of /i, and we 
look at their asymptotic behavior. 

The second means is a refined density mapping: we fix a measure defined by 
an initial density p{po] x). Then, we compute the Fourier coefficients Ck(y n ) via 
cq. (|2"Tj) . Of course, while p, n — > p weakly because of Mendivil's theorem, there 
is no guarantee that the densities p(p n ;x) will tend to a function. Therefore, 
we compute numerically the norms of p(p n ', x) reconstructed from their Fourier 
coefficients, and the distances in eq. (fT5|) . We then look for either divergence 
of the bounded variation norm (the L 1 norm must be conserved), or geometric 
convergence of the distances. In the first case we conclude for singularity of the 
measure p, in the second for absolute continuity. 

We can test this approach by examining three cases of I.F.S. measures with 
two maps, a — |(A_i + Ai), and S > \. The first two values of 5 are chosen 
in the two denumerable sets for which the spectral type is rigorously known: 
absolutely continuous in the first case, a) S = 2~ 1//2 ~ 0.707I0678II865 and 
singular continuous in the second case, b) 6 = l/pi ~ 0.7548776662467, p\ being 
a Pisot number. We also consider the case c) S = 3/4 = .75 which is pretty 
close to b). The last case is particularly interesting, since it is conjectured 
that for rational values of S larger than one half the measure y is absolutely 
continuous. In this specific case our analysis shows that the conjecture appears 
to be numerically validated. 

Observe Fig. [3] the analysis of the L 1 distances \\p(fi n ) — /?(/x n _i)||i versus 
n shows exponential decay (hence convergence in virtue of Cauchy criterion) for 
the known absolutely continuous case b) and also for the conjectured case c). 
In case a) one observes to the contrary (as expected) divergence of the bounded 
variation norm of p{\x n \x). Also notice that difference between cases b) and c) 
can only be appreciated after about ten iterations. 

We also plot in fig. [4] the absolute values of the Fourier coefficients Ck{p) 
of the three invariant measures, versus k, in double logarithmic scale. This 
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Figure 3: Bounded variation norms of densities p n for the three I.F.S. described 
in the text: a) open squares, magenta; b) filled squares, light blue; c) open 
circles, yellow. Also plotted are the L 1 distances between successive densities: 
a) green crosses; b) red pluses; c) blue asterisks. 



confirms the the results obtained by the analysis of the iterative algorithm. 
Coefficients for the case b) can also be computed analytically. 

7 Test Case III: a is a Singular Continuous Mea- 
sure 

We have seen so far two cases where a is a discrete measure or an absolutely 
continuous one. We want now to consider a singular continuous measure a. In 
order to be able to investigate this case numerically, we need a fast and reliable 
way to compute the Fourier transform of a. But indeed we have already at 
hand a measure a with these characteristics. Consider in fact the measure n 
generated, for 5 < 5, by the Bernoulli distribution, section[5]and[6] It is singular 
continuous and in addition its Fourier transform can be easily computed via 
Algorithm 2. Therefore, by taking now this measure to be the new distribution 
a, we can apply Algorithm 2, in a recursive fashion. 

The basic idea of this recursion is simple: suppose that the prey moves 
with distribution a((3) — /^°)(/J), this latter being an arbitrary measure. The 
predator is attracted to the prey according to the affine maps with contraction 
ratio Si: this determines its invariant distribution (jP-K In turn, this predator 
is hunted by a second species: for this last animal, a(/3) = fi^(f3) and 82 is the 
contraction ratio: its distribution is then /j( 2 \ Clearly, there is no limit to the 
number of species in this food chain: this fact has an interesting mathematical 
formulation that will be exploited elsewhere. For simplicity we investigate in 
this section only a simple case: //°) is the Bernoulli measure, 5% < |, and the 
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Figure 4: Fourier coefficients |cfc(/x)| vs. k for the three I.F.S. measures a) (red 
dots) b) (green dots) and c) (magenta dots) described in Sect. [6] Solid lines 
are drawn to mark the asymptotic behaviors fc -1 / 2 (borderline for assessing 
absolute continuity), k^ 1 and k~ 2 . 



choice cr(/3) = /j,( l '(j3) leads to an I.F.S. with uncountably many maps, whose 
fixed points populate a Cantor set. 

The theory developed so far enables us to examine numerically the invariant 
measures of these I.F.S. As the typical case, we observe that the measure /jffl 
is smoother than fx- 1 '. It can indeed be absolutely continuous even when jjS 1 ' 
is singular continuous and supported on a Cantor set. This case is pictured in 
Fig. [SJ in which //*) is generated by an I.F.S. with fj,^ — ^(A_i + Ai) and 
Si = 0.4, and we choose two values of the second contraction ratio: 62 = 0.1 
and 62 = 0.2. In both cases we observe absolute continuity of the measure ^ 2 \ 
with a continuous density. The densities of the two measures fit- 2 ) are plotted in 
Fig. [5j Observe that the support of the measure fi^ includes that of fi^ [9]. 
Indeed, a stricter "inequality" formula between them can be derived. We also 
plot in the inset the geometric convergence of the BV and L 1 distance functions 
(fT8|) . The initial increase of the bounded variation distance is due to the fact 
that we have employed here an initial gaussian distribution (see sect. [4j. 

Should one be tempted to conjecture that /i^ is always absolutely continu- 
ous in this setting, in the next and final section we shall meet a case when // 2 ) 
is singular continuous. 



8 A computer assisted proof of absolute conti- 
nuity 

In this final section, we show how a computer assisted proof of absolute conti- 
nuity of the measure fi can be developed in a specific case. 
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Figure 5: Densities p(/j,^;x) of the two I.F.S. measures described in the text. 
Labels are: 62 = 1 (green curve), 82 = -2 (red curve). For symmetry, only 
half of the curves are shown. In the inset, bounded variation and L 1 distances 
||p(// n +i) — p{Hn)\\ versus iteration number n: labels are 82 = .1, BV (green 
curve, crosses), L 1 (magenta, squares); 62 = -2, BV (red, pluses), L 1 (blue, 
asterisks) . 

Lemma 8.1 Let the measure /i^- ' with support in [— 1, 1] be given, and let fi^ 1 ' 
be the invariant measure generated by an affine, homogeneous I.F.S. with con- 
traction ratio 6 and with distributions of fixed points /i^ ' . Let also p\ 2 > be the 
invariant measure generated by a second affine, homogeneous I.F.S. with con- 
traction ratio 8 (the same as before) and with distributions of fixed points fj}- 1 '. 
Then, for any N > and for any t, the inequality holds: 



log \^ 2 ){S- N t)\ < N{\og \iiW(6t)\ + log \n^(S~ ^t)\} 



(23) 



Proof. Use cq. (|22|) to express \ogpt^(y), and use the convenient notation 
&j(y) '■= log \fi^(y)\. Gather terms (some algebra required) to obtain 



JV-l 



$ 2 (S- N t) = J2 ( l + ^a(S l - N 6 2 t) + mi(h) + * 2 (t). 



(24) 



1=0 



Observe that all functions $j are less than, or equal to, zero so that by keeping 
only the I — N — 1 term in the summation at r.h.s. of eq. ([24)) and discarding 
also the term $2(t) leads to the inequality (f23| . This is crucial for proving: 



Theorem 3 Let /J,^', for j — 0,1,2, be defined as in Lemma \8.1[ Let also 
r, = -sup-j^oOr 1 ^) + $i(<ft); t <= [l,^ 1 ]}. Then, if r) > ^log^ 1 ) the 
measure ir- 2 ' is absolutely continuous with an L 2 density, and if r\ > log((5 _1 ) 
the density p(/i; x) is continuous. 



13 



Proof. From Lemma 18.11 it follows that for all AT > 0, and for all t E 
[5~ N ,5~ N ~ 1 ] the function $2(*) is less than —r/N. One can therefore apply the 
Sobolev criterion described in Sect. [51 I 

This theorem quickly turns into a computer-assisted proof of absolute con- 
tinuity. In fact, notice that the supremum required to compute r\ is indeed the 
maximum of a continuous function defined on a finite interval. As such, it is 
easily computable with sufficient precision to establish whether the inequalities 
in Thm. [3] are satisfied. For the systems we have been considering in this paper 
this is easily accomplished, by computing the Fourier transforms via algorithm 
2. As just an instance of this procedure, we can prove 

Theorem 4 Let ij}^ = |(A_i+Ai), and let fj!^', j = 1,2 defined as in Lemma 
\8.1\ Then, if 5 = I the measure fj,^ 2 ' is absolutely continuous with a continuous 
density, and if S = jx it is absolutely continuous with an L 2 density. 

Proof. In the two cases quoted, one can easily compute that 77/ log(<5 -1 ) ~ 1.5035 
and T]/log(S~ 1 ) ~ .96838, respectively. I 

It appears numerically that in addition to being a rigorous inequality, for- 
mula (|23[) is a rather precise estimate of the asymptotic behavior of the Fourier 
transform of ijS 2 \ So, by letting now 8 = j^ we find that the value 77/ log(<5 _1 ) ~ 
0.2218 implies that the corresponding measure fi^ is singular continuous. This 
conclusion is also confirmed by the iterative analysis detailed in Sect. [5J 

9 Conclusions 

We have described numerical techniques to determine the continuity type of 
measures generated by affine, homogeneous I.F.S. with uncountably many maps. 
These techniques are inspired by the theory of dynamical systems. In turn, 
they suggest methods to obtain rigorous, and computer-assisted proofs: one 
should focus on the existence (or non-existence) of an attractive fixed point in 
a suitable space of densities, most conveniently under the bounded variation 
distance. Alternatively, one must consider the decay rate of Fourier coefficients, 
also obtained via fixed point procedures. 

We have sketched examples of both techniques, that have helped us to clarify 
the fact that such I.F.S. measures can be of all continuity types. According to 
the language adopted in this paper, we have described the implications of the 
continuity of a on those of fj, in what we believe to be a rather detailed pic- 
ture. We have also shown applications of our methods to the classical problem 
of infinite convolutions of Bernoulli measures, for which determining the conti- 
nuity type of a specific measure whose parameter does not belong to a known 
countable set of values is still an open problem. 
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